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Recently, some over-luminous la supernovaes are found, suggesting that their progenitors are white 
dwarfs more massive than the Chandrasekhar limit, which perhaps result from ultra-strong magnetic 
field inside the white dwarfs. We present an equation of state, explicitly magnetic-dependent and 
analytically practicable, and observe that the change of equation of state due to magnetic field 
waning along radius will so significantly influence the configuration of a white dwarf that its density 
does not monotonically decrease, but goes down at first, re-peaks near the crust and falls off again. 
As a supernovae will, in single degenerate la supernovae system, leave the remnant of its companion 
and a neutron star (pulsar star), we point out that the observations of these objects can put our 
model into tests. 


I. INTRODUCTION 

A main sequence star of low or mediate mass will fi¬ 
nally evolve into a white dwarf (WD) after supernova 
explosion. S. Chandrasekhar[1] has shown that a WD 
will not sustain a mass over 1.4 Mq, at which very mass 
a WD, after absorbing mass from its companion, will 
explode again and become the la Supernovae (SNIa). 
This is the Single Degenerate (SD) mechanism of SNIa 
formation|2J[3], and suggests that SNIas are all of with 
nearly equal luminosity and hence can be regarded as the 
cosmic standard candle to measure cosmology distance, 
contributing to the surprise discovery of the accelerative 
expansion of our observable cosmos [3], and revolutioniz¬ 
ing our understanding of fundamental physics world even 
radically extending our conventional gravity and cosmol¬ 
ogy models. However, some SNIas, such as SN 2007if, 
are recently found whose progenitors are even twice as 
massive as Chandrasekhar limit (ChL)[5]IS][7]IH]i[ini, 
for which some authors try to introduce electric field pT] 
or rotation contributions jT2( into a WD to bypass ChL 
within the SD SNIa model. Besides, Double Degener¬ 
ate (DD) Mechanism [I3J, where SNIa is ignited by two 
WDs merging together, can account for a wide range of 
SNIa luminosity while fails to explain why most SNIas 
are similar. 

In Chandrasekhar’s original work[T] , the equation of 
state (EOS) is determined by the degenerate pressure 
amongst relativistic Fermi gas of electrons. To accom¬ 
modate the recently found super-massive DWs, the work 
by U. Das and B. Mukhopadhyay (D&M) 10 !£>' issues 
that, since ultra-strong magnetic field (USM) will quan¬ 
tize free electron gas into Landau levels and soften the 
EOS, a WD will be 2.58Mq’s massive when inside is with 
a constant magnetic field as strong as 10 18 Gawss(G). 
This work is, however, challenged for over-simplicity and 
on stability problems [16] [17] [18]. Not very long, U. Das 
and his collaborators[Dj refine their previous workfH] by 
relativistically considering a density dependent isotropic 
magnetic field, and response some critics [20j}21j. Mean¬ 
while, R. Nityuananda and S. Konarj22j attempt to con¬ 
strain the configuration of a magnetic WD but their re¬ 


sult was soon proven erroneous }23j. 

The EOS used in the above papers is simply a 
quadratic approximation for only v = 1 Landau Level 
in USM existence, because the original EOS with Lan¬ 
dau levels is in the awkwardly complicated form of al¬ 
gebraic summation. As the magnetic field will decrease 
with radius, the EOS will surely shift accordingly, which 
is, however, overlooked. We will reduce the Landau level 
EOS to a feasibly analytical form without loss of accu¬ 
racy, which will simplify the discrete summation and ac¬ 
curately demonstrate the influence of varying magnetic 
field on the multi-level EOS[14| |24j. 

General Relativity contains high non-linearity from 
complication of curvature, with appearance of metric ten¬ 
sor in energy-momentum tensor, and even interlocked 
couple of electromagnetism in the magnetized astro¬ 
physics. Especially, the structure equation will be in¬ 
surmountable to solve exactly when magnetic field de¬ 
stroys the simple spherical symmetry. In D&M, they cir¬ 
cumvent such complexity by ideally considering constant 
magnetic field in Newtonian frame. Some authors also 
try to reduce the complexity by unnaturally introduc¬ 
ing into a spherically symmetric metric with anisotropy 
in energy momentum tensor[25], Only few works |26j [27] 
deal with non-spherical cases by virtue of, for example, 
GRMHD code|2Hj, but the magnetic fields in their works 
are much too weak. It is well-know that as the gravi¬ 
tation field inside a white dwarf is weak [29] , the linear 
theory will be sufficient enough for describing. Hence, we 
are going to develop the general equation of equilibrium 
in the framework of linearized General Relativity, where 
anisotropy property can be treated much easier. 

We arrange this paper as follows. The next section 
will present the linearly relativistic equation of equilib¬ 
rium for our working frame. The third section will be 
devoted to magnetic field distribution by modelling with 
two regions, and the forth section deals with the equa¬ 
tion of state in the strong magnetic field. Finally, we 
will summarize the result and draw the conclusion with 
some discussions for possible perspectives. Throughout 
this paper, the speed of light c = 1 is taken, and the 
Minkowski metric is set to be diag (—1,1,1,1) as com- 
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monly adopted. 


II. LINEARIZED GENERAL RELATIVITY 

Denote the deviation from Minkowski metric as, 

"Uiv = 9nv ~ ( 1 ) 

so that the Einstein’s field equation, expanded to the 
linear order and imposed the Lorentz gauge, will be 

□7 ilv = -167tG (t^ - . (2) 

If inserting into T MJ , as summed from that of ideal fluid 
and electromagnetic field, noticing that T = 3p — p, re¬ 
garding ^700 as the potential </>, we get the modified Pois¬ 
son Equation, 


A (f> = 47 tG (p + p B + 3p), (3) 

where pb = is the energy density of magnetic field. 

We ignore the rest part of 7 ^ for they will not appear 
in the equilibrium equation. 

The equation of motion can not be directly obtained 
from linear theory as = 0 ; however, for ideal fluid 

in the frame of General Relativity, if we assume that there 
is no 4-current inside a star, the equilibrium equation is 
given as 3[)| 


Vp = -(p + p)VlnV-3oo, (4) 

which, expanded to the leading order, will be 

Vp = ~{p + p)V<j>. (5) 


These equations will be formed as a close system if 
blessed with the equation of state and magnetic distri¬ 
bution. 

For illustrious purpose, consider the spherically sym¬ 
metric case, in which Eqs. ([3]) and ([ 5 ]) will be reduced 
as, 


d M 

dr 

dp 

dr 


4wr 2 (p + p B + 3p) 


~{P+P) 


GM 

r 2 


III. MAGNETIC FIELD 


( 6 ) 

( 7 ) 


initially illustrate the effect of a constant field, which, for 
magnetic field will dwindle at least two orders from cen¬ 
ter to surface, are only a simplified toy model with some 
debattings. 

However, a well-known magnetic field distribution 
is indeed constant in the center core while decreasing 
rapidly out of the core. As it is known that for an 
uniformly magnetized sphere in vacuum, magnetic field 
is constant inside the sphere and dipole distribution 
outside [34] , 


B = 


Bo 

— ^V(B„.r/r») 


r ^ r 0 
r > r 0 


( 8 ) 


where r 0 is the radius of inner core. 

Because the main purpose of this paper is to illustrate 
how the impact of decreasing magnetic field on EOS will 
affect the structure of an USMWD, we will assume that 
the magnetic field is spherically symmetric and, in light 
of the distribution of uniformly magnetic sphere, we take 
the magnetic field as constant inside an inner core in the 
center of a WD and inversely cubic outside, that is, by 
generalizing the above distribution directly 


b= B “ 313 r ^° ( 9 ) 

|B 0 "ro/ r r>r 0 

where ro is the radius of the inner core. The direction 
of magnetic field induction will accommodate to satisfy 
Maxwell’s equation. 

In addition, the conservation law of the magnetic flux 
also suggests that the ruminant neutron star of an SNIa 
possess a magnetic field of at least 10 2O G. As neutron 
stars are much easier to detect, it is very likely to dis¬ 
cover USM neutron star by the pulsar star properties 
like the timing. This will justify the USMWD hypothesis 
whether the neutron star is from a SNIa or not, because 
if it is, its progenitor’s magnetic field will be sufficiently 
strong; if not, at least the existence of a main sequence 
star of magnetic field stronger than 10 lo G will be con¬ 
firmed, which also can be the progenitor of an USMWD. 


IV. EQUATION OF STATES 

Landau level describes a free electron moves in an uni¬ 
form magnetic field B. Set the magnetic field along z 
direction, i.e., the vector potential be 

Afj = (0, By, 0,0); (10) 


The magnetic field inside an ultra-strongly magnetic 
WD (USMWD) is not clearly known. As the mag¬ 
netic flux BR 2 is preserved through a collapse process 
[3Tj [32] [33] . an usual DW can succeed from its progeni¬ 
tor a magnetic field at least 5 orders stronger. However, 
accretion might also play an indispensable role in the for¬ 
mation of an USMWD. Das and his collaborators[14] [IS] 


so the Landau energy level, given by Dirac equation, is 
E u , Ps = VpI + m e { 1 + 2 vB d ) (11) 

where v = j + | -f- a is the quantum number of Landau 
level, and Bp = Bhe/m 2 . So the maximium of v is 
u m = (tp — l)/2 B d where the ep is the Fermi’s level, 
and the EOS will be El 
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( 12 ) 

(13) 



be surmounted by the famous Eular’s formula 
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(a )S P and its deviation 
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where Bk stands for the Ath Bernoulli number, as 1/6, 
1/30, 1/42, etc. and £ some number between 0 and 1. 
For £f > 1 and Bf> A> 1, the Landau EOS can be (for 
details, see Appendix A) 


P = 


P = 


2 p e m H B D 
(2tt) 2 A| 

|^K + 2% m -l)S> m )] 


\j2B D v m + Q(v m - 1 )S p {y rn ) (16) 

(17) 


Vm appears in the above EOS as a continuous param¬ 
eter because it sufficiently approximates the original dis¬ 
crete summation as illustrated in FIG. |T]. 


(b)S p and its deviation 

V. RESULTS 


FIG. 1. Exactitude of S p and S p . The black lines represent 
S p and S p and the red lines stand for the accurate original 
discrete summation. We can see that they are very close and 
get closer and closer as increases. The relative errors are 
plotted in dashed gray lines. 


where p e means molecular weight per electron, A e Comp¬ 
ton wavelength of electron, m F mass of Hydrogen atom, 
cf Fermi energy, 6(z) is the step function that is equal 
to constant one except vanishing when z negative, and 


Given the EOS, Magnetic field distribution, and struc¬ 
ture equations, we obtain the Mass-Radii relations as 
plotted in FIG. [2j Masses will increase with respect to 
core radii and central strength of magnetic field, while 
Radii are almost determined by the ro and independent 
of B d . Despite the magnetic field, p, p , and I'm. are also 
nearly uniform in the inner core, as for v m inside, pCp, 
Pb p, for its derivative, 


y/2p2 bt-h GM 
2m e B 1 0 /2 r 2 


3.3 x 10~ u B 0 r. 


(18) 


V 0 ) 



In [z + yj\ + z 2 ) 


(14) 


The summation terms that appear in EOS cause awk¬ 
wardness and inconvenience for practical use, but it can 


However, outside the inner core, the derivative of Bp 
will also contribute to v ' m , together with the accumula¬ 
tion of mass, so that v m will increase significantly and p 
will begin dropping, i.e., p keeps constant inside whereas 
decreases outside the inner core, and the later p dwin¬ 
dles, the slower p vanishes; hence the increasing of radii 
and mass with respect to r 0 . The influence of magnetic 
field strength can be understood as follows. M, as is 
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3 1‘l 

dominated by p, is proportional to at most B^ while 
p ~ Bp', therefore, p will increase faster than its deriva¬ 


tive. 


We will study a typical USMWD with Bd = 20000, 
Ro = 1000m, whose mass is solved to be 2.40Mq and 
radius 51.930 km. Its density and Vm distribution is pre¬ 
sented as in FIG. [3] . FIG. 3(a) indicates that p does 


not monotonously decrease with respect to r. In the con¬ 
trary it will peak again near the crust surface, which 
violates P. Bara’s negative result on the magnetic influ¬ 
ence structure [3B]. When v m increases as Bp goes down, 
p ~ p 4 / 3 B]f 3 . Accordingly, p will increase if p decreases 
slower than 1/r. Such structure undermines the Global 
stability as mentioned in Ref. mi , 


E, 


mag 


fo M fdm 


E, 


grav I Jo 1 


(19) 


In this expression, E mag « f Airr 2 psdr is nearly indepen¬ 
dent of p(r) whereas E grav = f M(r)M'(r)/rdr peaks at 
MqS(t) and wanes when p is not centralized near the 
center. The global instability for this DW is as high as 
60.2. 

In addition, as magnetic field generally introduces 
dipole effect, the dipole momentum of an USMWD will 
be amplified. If an USMWD is found in a binary system, 
the ” re-peaking” effect can be put to test by comparing 
the dipole mass momentum measured from the shift of 
the binary system period with both its mass and radii. 



Radium m 

(a )p — r relation 



Radium m 

(b)i 'm — r relation 


FIG. 3. Structure of a WF with Bd = 20000, and R0 = 
1000m 


VI. CONCLUSION 

We have derived the analytical form of Landau level 
EOS, presented the M-R relation for USMWD, and have 



Radii m 


FIG. 2. Mass-Radii Relations. The M-R relation of 
B d = 30000, 20000,15000,12000,10000, 8000, 6000 is drawn 
in brown, red, orange, gray, green, blue and purple. The 
dashed black lines in the figure represent the M-R relation for 
same inner core radii, from 500m to 2000m at the interval of 
250m. 


found the ” re-peaking” effect of the density profile for 
the EOS becoming harder when magnetic field wanes. 
We also prove that USMWD’s high magnetic field hy¬ 
pothesis suggests that the existence of highly magnetic 
neutron star and of SNIa ruminant companion. If they 
are both discovered, USM will be possible and SD SNIa 
model be confirmed. At least, we can conclude that the 
progenitors of some over-luminous SNIas will indeed be 
USMWD. We should say that the concrete mechanism 
and distribution of such USM are still unclear, which are 
however, indispensable for deriving the statistical distri¬ 
bution of USMWDs. This is only the first step for our 
work, though only touching some aspects of the inter¬ 
esting star forming processes involved and with upcom¬ 
ing more accurate astrophysics observational data we are 
keeping on studying the accretion process in order to un¬ 
derstand further the origin of the USM, as it may shed 
new light on our novel understanding the extrema condi¬ 
tional physics. 
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Appendix A: S p and S p 


In this appendix we will present more details for the 
analytical treatment to the complicated EOS. We take 
Eular’s formula to the order n= 1 and define E^ as, 

E E (F(i),l + l,m-l) 

r n i i 

= J F(x) dx - - [F(l) + F(m)} + - [F\m) - F'(V 

(Al) 


So S, 


pi 


s; « o + 1 + E1, - 2) 

_ -1 - 5y/Frn ~ 12 I'm + 16 ^ 

24 yJVm 


(A2) 


S p = S* + J 2 to secure that S p vanishes at = 1. 


Similarly, for S. 


pi 


pl(i) = 2v m - 


i - 


i — 2ty, 


—- -In Av m ( A4) 

vL 2 


p2(i) = 2 ln * 


(AS) 


s ;=E‘ 


-i+ 



- ln(w— 1 + 


i=l 


(pl(i) + p2(i)) 


~ E s (pl, l,i/ m - 1 ) +pl(u m ) 

+ E£)(p2, 2, v m — 1) + p 2(l) + p2(i/ m ) _ 

= \-Mu m + (315 - 372 In 2)i& -7 In V -f 
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